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The pairing symmetry is examined in highly electron-doped Ba(Fei_ccCoccAs )2 and AyFe 2 Se 2 (with A=K, 

Cs) compounds, with similar crystallographic and electronic band structures. Starting from a phenomenological 
two-orbital model, we consider nearest-neighbor and next-nearest-neighbor intraorbital pairing interactions on 
the Fe square lattice. In this model, we find a unified description of the evolution from s±-wave pairing (2.0 < 
n < 2.4) to d-wave pairing (2.4 < n < 2.5) as a function of electron filling. In the crossover region a novel 
time-reversal symmetry breaking state with s± id pairing symmetry emerges. This minimal model offers an 
overall picture of the evolution of superconductivity with electron doping for both s±-wave [Ba(Fei_a;Cox As) 2 ] 
and d-wave [AyFe 2 Se 2 ] pairing, as long as the dopants only play the role of a charge reservoir. However, the 
situation is more complicated for Ba(Fei_xCoccAs) 2 . A real-space study further shows that when the impurity 
scattering effects of Co dopants are taken into account, the superconductivity is completely suppressed for 
n > 2.4. This preempts any observation of d-wave pairing in this compound, in contrast to A^Fe 2 Se 2 . 

PACS numbers: 74.70.Xa, 74.20.Rp, 74.25.Dw 


The recent discovery of Fe-pnictide based superconductors 
offers a new family of materials in which the nature of super¬ 
conductivity can be explored ESI- Although the mechanism 
of superconductivity (SC) in this family still remains an open 
subject, several theoretical models indicate that the complex 
geometry of the Fermi surface (FS), including both hole and 
electron pockets in the Brillouin zone (BZ), should be mainly 
responsible for the SC with 5±-wave pairing symmetry in 
weakly and moderately doped systems IbUTOl. In the FeAs- 
122 family, the SC in the electron-doped Ba(Fei_a;COa;As )2 
disappears almost at the same doping level where the hole-FS 
pockets vanish near the F-point El. This feature has been 
widely interpreted to indicate the existence of a correlation 
between the SC and the FS topology of the Fe-pnictide com¬ 
pounds. On the other hand, a new series of iron-chalcogenide 
122-compounds such as A^Fe 2 Se 2 (A=K, Cs), have been dis¬ 
covered recently with relatively high SC transition tempera- 
ture (Tc) of about 31 K EHEl- These superconductors are 
the most heavily electron-doped among the iron-based com¬ 
pounds with 0.8 < ^ < 1. Electronic band structure calcula¬ 
tions for A^Fe 2 Se 2 indicate that only electron pock¬ 

ets exist near the M-point of the BZ. A recent angle-resolved 
photoemission spectroscopy (ARPES) experiment El also 
showed the presence of electron pockets around the M point 
and the near absence of a hole pocket around the F-point 
at ^ = 0.8. For weakly and moderately electron-doped 
Ba(Fei_a^Coa^As )2 with x < 0.4, the FS consists of both 
electron pockets near M points and hole pockets near the F 
point. These compounds are thought to have s±-wave pairing 
symmetry due to magnetic fluctuations between the M and F 
points in on the FS. However, for the heavily electron-doped 
A^Fe 2 Se 2 with ^ = 0.8 to 1, the FS contains no hole pockets 
at the F point and the proper superconducting pairing symme¬ 


try should not be s±-wave. Instead, it was predicted to have d- 
wave pairing symmetry (22^261. Therefore, it is very impor¬ 
tant to understand the evolution of the superconducting pair¬ 
ing symmetry within a real-space formulation for A^Fe 2 Se 2 
with 0.8 < ^ < 1 and at the same time why SC completely 
disappears for Ba(Fei_a^C 0 a:As )2 with x > 0.4. 

So far the s±-wave pairing symmetry in Ba(Fei_a^C 0 a:As )2 
has been attributed to a next-nearest-neighbor (NNN) pairing 
interaction among the electrons of the Fe atoms |[9l|23. In 
addition, its phase diagram has been mapped out as a func¬ 
tion of electron doping ||28H3Q1 and it has been shown that 
SC vanishes for electron doping x > 0.4. This interpreta¬ 
tion seems plausible, since at this point the hole pockets at 
the F point vanish as the hole states are completely fllled by 
doped electrons. Not only are these results consistent with 
experiments |[3Tl - [34ll . they also suggest that the NNN-pairing 
interaction is able to capture the essential ingredients of the 
magnetic fluctuations between the M and F points El. Ac¬ 
cording to the above theoretical works ll28H30l . there is no SC 
in A^Fe 2 Se 2 with 0.8 < ^ < 1, in clear contradiction with 
recent experiments. 

In this work, we seek to obtain a unifled picture of the evo¬ 
lution of the superconducting pairing symmetry with mod¬ 
erate pairing strength for the above mentioned compounds 
by including the nearest-neighbor (NN)-pairing interaction 
in addition to the NNN-pairing interaction I^nn- Since 
the electronic structure of A^Fe 2 Se 2 is similar to that of 
Ba(Fei_a^C 0 a^As) 2 , we will employ this new Hamiltonian to 
re-examine the SC of both materials. The strength of the 
pairing interactions is chosen according to two factors, one 
is that the pairing symmetry of Ba(Fei_a^Coa;As )2 is domi¬ 
nantly 5±-wave at doping level x < 0.4, the other one is 
that the A^Fe 2 Se 2 has a flnite SC order parameter at dop- 
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ing level 0.8 < y < 1. The Hamiltonian is numerically 
solved on a 28 X 28 lattice by using the Bogoliubov-de-Gennes 
(BdG) equations. Our results reveal that the pairing symme¬ 
try changes from 5±-wave to d-wave as the electron doping is 
increased. For moderately electron-doped Rb^Fe 2 Se 2 (with 
y = 0.3), the experimentally observed SC |[35]| should have 
an s±-wave pairing symmetry. In Ba(Fei_a^COa^As )2 com¬ 
pound, although the SC is predicted for x > 0.4 in the present 
calculation, it has not been observed in experiments MM- 
We will numerically demonstrate that the d-wave SC in this 
region can be completely suppressed by scattering due to ran¬ 
domly distributed Co atoms in the Fe planes. 

We use a phenomenological two-orbital tight-binding 
model tSOl to numerically perform calculations of the FS, lo¬ 
cal density of states (LDOS), magnetic and SC order parame¬ 
ters. This model has been successfully used to theoretically 
describe the generic phase diagram and other properties of 
FeAs-122 superconductors 1^ 13614381 . Based on the fact that 
they share similar crystal and band structures with the FeAs- 
122 family, it is reasonable to apply our two-orbital model, 
which is built to describe the FeAs-122 family, to A^Fe 2 Se 2 
compounds, after an additional NN-pairing interaction is con¬ 
sidered. 

Consider the Hamiltonian H = Hq ^ Hsc ^ Hint that 
describes the energy of charge carriers. Without the impurity 
term, here H follows the same formulation as Ref.|30|for the 
hopping (Ffo), pairing (Hsc) and on-site interaction (Hint) 
terms. We express the matrix form of H with the basis, ^^i^ = 
^ to calculate the 

eigenvalue and eigenvectors of 
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here both NN and NNN intra-orbial pairing orders are calcu¬ 
lated from the following equations |[30l . 
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where e G {NN, NNN} to denote the nearest-neighbor and 
next-nearest-neighbor pairing bond. Throughout this work, 
we use the same hopping and on-site interaction as in Ref.[3Ql 
The chemical potential is determined self-consistently cal¬ 
culating the averaged electron filling n from the BdG equa¬ 
tions. is the intra-orbital pairing interaction between NN 
sites, and Vnnn is the intra-orbital pairing interaction between 
NNN sites. All energies are measured in units of nearest- 
neighbor intra-orbital hopping jtsj, see definition in Ref. [3Ql 
The collinear spin density wave (SDW) order parameter is de- 
fined as = (-1)*" | “ (nini)), where the 

is the lattice number along x-axis and (n^^cr) is the electron 
density for site-i orbital-// and spin-cr. The pairing order pa¬ 
rameters of NN d-wave and NNN 5±-wave are defined as 

• = 8^1 J and 

• AE = 8^ respectively, 

where j = i ± is the NN sites of site i and = i±x ±y 
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FIG. 1. (Color online) SC and AFM phase diagram at T = OK for 
uniformly doped AyFQ 2 SQ 2 and Ba(Fei-ccCoxAs )2 as a function of 
the averaged electron number per site n, with NN pairing IAn inter¬ 
action chosen to be (a) 1.1 (b) 1.2 and (c) 1.3. Black squares, red and 
blue triangles represent the collinear SDW order parameter, pairing 
order parameter of NNN s±-wave and NN d-wave, correspondingly. 
The NNN pairing interaction is set to be fixed at Vnnn = 1.05. The 
purple shaded region represents the time-reversal breaking pairing 
state, s± id. 


is the NNN sites of site i, ex(y) = — i) • x(y) and N is the 

number of Fe lattice sites. We can write down the correspond 
form factor in k-space for and A'^^ under 2-Fe per unit 
cell Brillouin zone 

• Ad(k) = 2A^[sin(/ca;) sin(/c^)] and 

• A'^^ (k) = 2A'^^ [cos(k:jo) + cos(ky)]. 

In two of our previous works 1291 [361 . it was shown that 
the s±-wave superconductivity in the phase diagram for 
Ba(Fei_a^C 0 a;As )2 can be constructed from the NNN intraor¬ 
bital pairing interaction. The approach in these works 1^1^ 
gives a unified description only for s±-wave pairing of the en¬ 
tire phase diagram covering both the electron and hole doped 
regimes. The k-space s±-wave SC order parameter vanishes 
in highly doped regime (n > 2.4) where the hole pocket on 
the FS at the F-point shrinks to zero. On the other hand, SC 
has been observed in recently discovered Aa;Fe 2 Se 2 at higher¬ 
doping level than n = 2.4 1T214T51 with a predicted d-wave 
pairing symmetry, has been studied under k-space band pic¬ 
ture An insight of real-space picture of the d-wave 

pairing in addition to what was used in the previous works 
should be taken into consideration. 

In order to seek a unified theme in the real-space lat- 
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tice model for the pairing symmetry of the electron doped 
A^Fe 2 Se 2 and other FeAs-122 compounds, we find the NN 
intraorbital pairing interaction is the last piece of jigsaw to 
complete the entire picture. The phase diagrams thus obtained 
are shown in Fig. for fixed I^nn and several The av¬ 
erage electron-number per site is defined as n = 2 + x = 
2 + ^/2. The coexistence of SC order parameters from differ¬ 
ent pairing symmetry are labeled in the figure, together with 
2x1 collinear SDW order. In the calculation the NNN pair¬ 
ing is set to be constant Vnnn = 1-05, but the NN pairing 
is changing from Vnn = 1.1 to 1.3. When Vnn = 1.1 as 
shown in Fig. [Ja), from n = 2.0 to n = 2.4, the pairing 
symmetry is dominantly s±-wave and the d-wave component 
is very weak, while in the doping region 2.4 < n < 2.5, the 
SC has a dominatly d-wave like symmetry originating in the 
NN pairing interaction. There exists a sharp transition of the 
SC pairing symmetry at n = 2.4. When Vnn is increased 
to 1.2 as shown in Fig.[2b), the d-wave dominant region ex¬ 
pands down to n = 2.38 where pairing symmetry becomes a 
mixed complex one, s± -h id (2.38 < n < 2.4), meanwhile 
the magnitude the pure d-wave SC order increased a bit in 
the highly electron doped region (2.4 < n < 2.5). The SC 
and AFM orders remain unchanged at all other doping levels. 
When the NN pairing interaction is further increased to 1.3 as 
shown in Fig.[2c), the complex d-wave SC order further ex¬ 
pands to lower doping at n = 2.36, but still enhanced in the 
highly electron-doped region. The magnetic order and the SC 
pairing symmetry in Fig. [Jc) at other doping levels remain 
the same as the previous two figures. In both the Vm = 1.2 
and 1.3 cases, the regions for pairing symmetry transitions be¬ 
come broadened and the pairing symmetry becomes a mixed 
complex one, s± -f id, for 2.36 < n < 2.4. If the NN pairing 
interaction is larger than 1.3, the d-wave pairing order will 
emerge in optimal-doped region or even in the under-doped 
region. This is clearly not in agreement with experiments per¬ 
formed on Ba(Fei_a^Coa;As) 2 . It is also not reasonable to 
choose < ^nnn = 1-05, because the obtained d-wave 
SC order parameter in A^Fe 2 Se 2 (with ^ ^ 0.8 to 1) would 
be too small to explain the experiments. 

The FS topology is also responsible for the SC pairing 
symmetry of the system. The zero-temperature FS is de¬ 
fined by zero energy contours of the quasiparticles, which 
can be drawn by using the Fourier transformation of the mini¬ 
mal hopping Hamiltonian. In highly electron-doped samples, 
the SDW order is completely suppressed and thus we show 
the corresponding FSs in the BZ with two Fe atoms per unit 
cell. In the following calculation we choose = 1-1 and 
= 1-05, the same parameter as shown in the phase di¬ 
agram Fig.[2a). Fig.j^a) shows the FS in optimal-doped re¬ 
gion, where n = 2.15. Here the hole pockets around F-point 
located at the center of the BZ can be clearly seen. This result 
is very similar to that in a previous work ||36l, which stud¬ 
ied the evolution of the FS topology in FeAs-122 compound. 
Fig.j^b) shows the FS at n = 2.4, where the s±-wave SC is 
suppressed to zero and the hole pocket on the FS shrinks to 
a point at the center of the BZ. We find that there is no hole 



FIG. 2. (Color online) The band structure and its corresponding 
Fermi surfaces at doping level n = (a) 2.15, (b) 2.4, and (c) 2.46 in 
the Brillouin zone of with two Fe atoms per unit cell. Black(red) lines 
represent the electron(hole) Fermi surface, and the orange (dashed) 
lines in (a)/(c) represent the nodal lines of s± / d-wave SC order 
^'s± {k) ! A(i(/c) which are not crossing the Fermi surface. The plus 
and minus sign indicate the sign of the SC pairing order parameter in 
each region. 


pocket structure for all doping levels higher than n = 2.42. 
Fig. [^c) shows the FS at n = 2.46, where k-space d-wave 
pairing order parameter Ad{k) = 2Adsm{kx) sm{ky) is 
maximized. The orange dashed lines here denote the nodes 
(or zeros) of Ad{k) and they do not cross the FS. In another 
word, Ad{k) is positive on the FS near (tt, tt) and (—tt, —tt), 
and negative near (— tt, tt) and (tt, —tt). The sign change over 
neighboring electron pockets demonstrates that the SC phase 
from n = 2.4 to 2.5 shown in Fig.[2a) is of d-wave symmetry 
without nodes. Our FS calculation result is in good agreement 
with a random phase approximation (RPA) calculation 1^ . 
It should be noticed that the pairing symmetry changes from 
5±-wave to d-wave at the same doping level where the hole 
pocket disappears on the FS, which reveals the intrinsic cor¬ 
relation between FS structure and system’s pairing symmetry. 

We now proceed to study the effect of real-space inhomo¬ 
geneity and address the outstanding question as to why the 
SC has not been observed in Ba(Fei_a^COa^As )2 for n > 2.4 
where the SC has a d-wave symmetry according to the present 
calculation. In this compound, there should be in addition to 
charge doping also significant scattering of the itinerant elec¬ 
trons due to randomly distributed Co atoms in the FeAs layer. 
The disorder concentration can become rather densed when 
the sample is in the highly (Co) doped region. We specu¬ 
late that the densed disorder scattering suppress the d-wave 
SC in this compound. On the other hand, the Fe-planes in 
A^Fe 2 Se 2 are quite clean because the doped A-atoms are be¬ 
tween FeSe-layers, thus the impurity potential of the A atoms 
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FIG. 3. (Color online) The phase diagram at T = OK of 
Ba(Fei-xCoxAs )2 as a function of n after averaged over the ran¬ 
domly distributed Co impurities. Black squares, red and blue tri¬ 
angles represent the collinear SDW order parameter, pairing order 
parameter of NNN s±-wave and NN d-wave, correspondingly. The 
impurity strength is set to be Vimp = —1, the NN and NNN pairings 
are set to be Vnn = 1.1 and = 1.05. The gray shaded back¬ 
ground represents the boundary of the impurity-free calculated phase 
diagram which is the same as Fig.[2a). 

has little effect on electrons in the Fe-planes. This is why 
the (i-wave SC survives in highly electron-doped A^Fe 2 Se 2 . 
To consider the scattering effect of the disordered Co impuri¬ 
ties in Ba(Fei_a^C 0 a;As) 2 , the Hamiltonian due to the impurity 
part can be written as 

Himp — per-I (3) 

I pa 

where Vimp is the impurity strength at the J-th Co site in the 
lattice, the summation is over all randomly distributed impu¬ 
rity atoms. In this work, the impurity potential of Co is known 
to be weaker than those of Ni or Cu 1401 , and we set the im¬ 
purity strength to be Vimp = ~1- 

Choosing the NN and NNN pairing to be = 1-1 and 
F^nn = 1-05, the same parameter used to calculate the phase 
diagram in Fig.[2a), we have calculated the SDW and SC or¬ 
der parameters as a function of n by averaging over 10 differ¬ 
ent impurity configurations on a 28 x 28 lattice. Our results for 
the phase diagram of Ba(Fei_a^Coa^As) 2 , after averaged over 
10 impurity distribution configurations, are shown in Fig.as 
a function of n = 2 + x with x as the concentration of doped 
electrons or Co impurities. Since the statistical error bar at 
each point is smaller than the symbol itself, they are not shown 
on the graph. Although the SC with 5±-wave pairing symme¬ 
try still exists in the region for n < 2.34, the SC in the highly 
electron-doped region (n > 2.34) is completely suppressed 
by the disordered Co atoms. Here the (i-wave SC exhibited in 
Fig. for n > 2.4 is destroyed and Andreev bound states are 
created by the impurities. Similar phase diagrams are also ob¬ 
tained for Fnn = 1.2 and 1.3, but we do not show them here. 
The essential feature shown in Fig.[^is consistent with exper¬ 
iments on Ba(Fei_a^COa;As) 2 . The predicted (i-wave pairing 
symmetry for A^Fe 2 Se 2 with 0.8 < ^ < 1 is very robust even 
when F^n > 1-3. 


It is important to point out that the crystal structure of 
A^Fe 2 Se 2 with ^ = 1 is identical to that of Ba(FeAs) 2 . With 
1 > ^ > 0.8, minor disorder is introduced into the A-layers. 
But for ^ = 0.3, the crystallographic structure may greatly 
deviate from that of Ba(FeAs )2 and may not have a well- 
defined stable structure. Since the FeSe layer is not affected 
by A-atom doping, the SC observed for ^ = 0.3 should 
have s±-pairing symmetry. We predict that in A^Fe 2 Se 2 with 
^ = 0.8 to 1, the pairing symmetry should be dominantly d- 
wave and the SC is of considerable strength if Fnn > ^nn- 

In this work, we presented a unified description of the evo¬ 
lution of superconductivity by including both the NN and 
the NNN intraorbital pairing interactions. We showed that 
by starting with a phenomenological two-orbital model the 
pairing symmetry transforms from s±-wave to (i-wave in 
highly electron-doped A^Fe 2 Se 2 as its Fermi surface topol¬ 
ogy changed. The transition occurred when the hole pockets 
vanished near the F-point of the Brillouin zone as the hole 
states were completely filled by doped electrons. We also 
found the emergence of a complex pairing state, s± + id, in 
our calculations. The existence of such a time-reversal sym¬ 
metry breaking s±^id pairing state was suggested in previous 
studies from a pure band picture in k-space ||42l|43l. However, 
we need to point out that our real-space formalism allows the 
study of the local spectra, as well as the stability and robust¬ 
ness of these competing phases in the presence of impurity 
scattering states, e.g., at very low impurity concentrations, or 
of vortex core states in magnetic fields. We attribute the ab¬ 
sence of (i-wave superconductivity in highly electron-doped 
Ba(Fei_a^COa;As )2 to the presence of randomly distributed Co 
impurities in the Fe square lattice with weak scattering po¬ 
tential. In that sense Co doping is less destructive to super¬ 
conductivity than Zn doping, which suppresses the s± pairing 
symmetry at around 10% Zn concentration |[30l . 

Besides the overall evolution of the superconducting sym¬ 
metry from s±-wave to (i-wave pairing with doping, the other 
key result of our work is the detrimental effect of weak impu¬ 
rity scattering potential on superconductivity, namely at high 
Co concentrations. Until now the Co-impurity scattering ef¬ 
fect was ignored in calculations, mostly due to the nature of 
its weakness. Here, we demonstrated that weak impurity scat¬ 
tering may play a crucial role for explaining the absence of 
(i-wave pairing in some of the 122 iron-based superconduc¬ 
tors at high electron doping. Finally, we suggest experiments 
to probe the signature of the time-reversal symmetry breaking 
s± -h id state in highly electron-doped samples, however, with 
impurity-free Fe planes. The most promising region in the 
phase diagram of our two-orbital model is therefore around 
electron filling of 2.36 < n < 2.4, where small hole pockets 
exist. 
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